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In this study, we investigate the electronic and magnetic properties of graphane nanoribbons. 
We find that zigzag and armchair graphane nanoribbons with H-passivated edges are nonmagnetic 
semiconductors. While bare armchair nanoribbons are also nonmagnetic, adjacent dangling bonds 
of bare zigzag nanoribbons have antiferromagnetic ordering at the same edge. Band gaps of the H- 
passivated zigzag and armchair nanoribbons exponentially depend on their width. Detailed analysis 
of adsorption of C, O, Si, Ti, V, Fe, Ge and Pt atoms on the graphane ribbon surface reveal 
that functionalization of graphane nanoribbons is possible via these adatoms. It is found that C, 
O, V and Pt atoms have tendency to replace H atoms of graphane. We showed that significant 
spin polarizations in graphane can be achieved through creation of domains of H-vacancies and 
CH-di vacancies. 

PACS numbers: 73.22.Pr, 75.75.-c, 61.48.-c, 62.23.Kn 



I. INTRODUCTION 

After the synthesis of two dimensional (2D) 
graphene[T]-[5], its nanoribbons (NRs) have been a 
subject of interest, both experimentally [7, 8 and 
theoretically. [9UI5] Advances in experimental techniques 
are paving the way for integrating the exceptional 
electrical, optical and magnetic functionalities of 
these nanometer-sized materials into future electronic 
technology. In addition to numerous theoretical and 
experimental works on graphene and its NRs, research 
efforts have been also devoted to synthesize its various 
derivatives. Notably, the synthesis of a 2D hydrocarbon 
in honeycomb structure, [16 namely graphane, followed 
its prediction through theoretical works. [TTl ITS] Detailed 
analysis of hydrogenation processes of graphene leading 
to graphane and the existence of hydrogen (H) frustra- 
tions were discussed very recently. [15] Furthermore, it 
was reported that graphane NRs have more favorable for- 
mation energies than experimentally available graphene 
ribbons. [20J Recently, we reported the possibility of 
obtaining magnetization through dehydrogenation of 
domains on 2D graphane. [5T] Modification of electronic 
structure of graphane by introducing either a hydroxyl 
group or a H-vacancy was also investigated with GW 
self-energy calculations. [22] Although there are several 
recent studies on graphane, electrical and magnetic 
properties of its NRs have remained unexplored. 

Recent advances in graphane have motivated us to ex- 
plore the properties of zigzag and armchair graphane 
NRs. In this paper, using first-principles plane wave 
method within the density functional theory (DFT) we 
investigated the electronic and magnetic properties of 
bare and H- passivated graphane NRs. We also explored 
the effects of specific imperfections on these properties. 
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These are various foreign atoms adsorbed on the sur- 
faces of graphane NRs, vacancies and edge roughness. We 
found that these imperfections can attribute interesting 
functionalities by modifying the electronic and magnetic 
properties of graphane NRs. 



II. METHODS 

We perform first-principles calculations |23_ 
within DFT using projector augmented wave 
(PAW) potentials [24 and approximate exchange- 
correlation functional by spin polarized local den- 
sity approximation[25 (LDA). Kinetic energy cutoff 
h 2 \\t + G\ 2 /2m for plane-wave basis set is taken as 500 
eV. In the self-consistent potential and total energy 
calculations of 2D graphane a set of (35x35x1) k-point 
sampling is used for Brillouin zone (BZ) integration. 
The convergence criterion of self consistent calculations 
for ionic relaxations is 10~ 5 eV between two consecutive 
steps. By using the conjugate gradient method, all 
atomic positions and unit cells are optimized until the 
atomic forces are less than 0.03 eV/A. Pressures on the 
lattice unit cell are decreased to values less than 1.0 
kBar. 

In order to correct the energy bands and band gap val- 
ues obtained by LDA, frequency-dependent GWo calcu- 
lations are carried out. [26] Screened Coulomb potential 
W, is kept fixed at initial DFT value Wo and Green's 
function G, is iterated five times. GWo self-energy cal- 
culations are carried out with 15 A vacuum spacing, 
default cut-off potential, 160 bands and 64 grid points. 

Graphane NRs are treated by the supercell geometry 
within periodic boundary conditions. Vacuum spacing of 
at least 15 A is placed between adjacent graphane NRs 
to hinder the interactions. In specific cases double unit 
cell is used in our calculations to allow the possible anti- 
ferromagnetic (AFM) orderings along the ribbon edges. 
Reciprocal space integrations are carried out with 1x1x15 
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FIG. 1: (Color online) (a) Top and perspective view of the 
atomic structure of infinite 2D graphane sheet having honey- 
comb structure. Two sublattices of graphane are indicated by 
A and B. Black (dark) and blue (light) balls are for carbon 
and hydrogen atoms, respectively, (b) Charge density con- 
tour plots of diamond and graphane are shown on a plane 
passing through C-C-C-C and H-C-C-H bonds, respectively. 
The tetrahedral angle of the diamond 6c = 109.47°. Arrows 
indicate the direction of increasing charge density. The cal- 
culated values of 9c and 9h, namely C-C-C and H-C-C bond 
angles in graphane respectively, are given in Table [I] Contour 
spacings are 0.0286 e/A 3 . (c) The LDA energy band struc- 
ture where the orbital character of specific bands is also given. 
The band gap is shaded yellow/gray. 



Monkhorst-Pack k-point grids. 



III. TWO DIMENSIONAL GRAPHANE 

For a better understanding of graphane NRs we start 
with 2D graphane which is derived from the hydrogena- 
tion of graphene, where each carbon atom is saturated by 



a single hydrogen atom. Accordingly, the primitive cell of 
graphane consists of two carbon and two hydrogen atoms. 
Chair like configuration of infinite sheet of 2D graphane 
is formed, whereby each carbon atoms of A- and B-type 
sublattices are saturated by hydrogen atoms from above 
and below, respectively, as described in Fig. [ija). This 
is known as the most favorable and stable hydrocarbon 
conformation. [T71 I21j The planar graphene honeycomb 
structure, which is stabilized by planar sp 2 -hybrid or- 
bitals and 7r-bonding between adjacent perpendicular p z - 
orbitals is puckered (buckled) as a result of the adsorp- 
tion of H atoms, whereby a single C atom is bound to 
three adjacent C atoms and a single FT atom through 
'tetrahedrally coordinated' sp 3 -like hybrid orbitals. The 
buckling, 5, i.e. the perpendicular distance between A- 
type and B-type carbon sublattices, is calculated to be 
0.45 A. As a result, the lattice constant of graphene in- 
creases from 2.46 A to 2.51 A and hence the C-C and 
C-FI bond lengths become 1.52 and 1.12 A, respectively, 
as indicated in Fig. Ilja). The charge density contour 
plots given in Fig. [TJb) indicate that the charge distribu- 
tion in graphane especially around the C-C bond is remi- 
niscent of that in tetrahedrally coordinated diamond. In 
graphane, the angle 6jj between FT-C and C-C bonds and 
9c between adjacent C-C bonds are 107.35° and 111.51°, 
respectively. The mean value of these angles is equal to 
the tetrahedral angle 9c of diamond. The maximum of 
the C-C bond charge in diamond and in graphane are 
found to be ~ 0.282 and ~ 0.286 e/A 3 , respectively. 

The cohesive energy of graphane (per unit cell) relative 
to free C and FT atoms is obtained from E c = 2 x +2 x 
E<f - E^ ra , where E% and E$ are the total energies of 
single free C, and free FT, whereas E!f ra is the total energy 
of graphane. The cohesive energy depends whether the 
energies of magnetic or nonmagnetic (NM) states of free 
C and free H atoms are considered. Here the cohesive 
energy per unit cell of graphane is calculated to be 23.57 
(27.65) eV by considering the magnetic (nonmagnetic) 
states of free atoms. As for the average C-FT bond energy, 
one can use the formula E C -h = {-E<f ra + E% re + 2 x 
E^ )/2 where Er£ re is the total energy of graphene. The 
calculated value is 2.84 (3.74 for nonmagnetic case). We 
reported|21j that the desorption of a single FT atom from 
the graphane is an endothermic reaction with an energy 
of 4.79 eV. When compared with the energy of C-FT bond, 
the desorption energy of single FT is larger, since the single 
FT removal creates reconstruction of nearby atoms which 
further reduces the total energy of the system. 

Due to sp 3 -saturation of C atoms, 2D graphane is a 
NM semiconductor with a direct band gap of 3.42 eV 
as shown in Fig. [TJc). However, this band gap, which 
is underestimated within DFT, is corrected by GWo ap- 
proximation to become 5.97 eV. |21l |2"2"] While the top of 
the valence band originates from p x - and pj,-orbitals of C 
atoms, the bottom of the conduction band has C-p z or- 
bital character. Earlier, we also calculate |2T| the phonon 
modes of infinite 2D graphane, which yield all real fre- 
quencies in BZ. Having all frequencies positive indicate 
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TABLE I: Comparison of the calculated quantities of graphene and graphane. Lattice constant, a; C-C bond distance, dc—C, 
C-H bond distance, dc-H', the buckling, S [see Fig. [TJa)]; angle between adjacent C-C bonds, 6c [see Fig. [TJb)]; angle between 
adjacent C-H and C-C bonds, Oh; energy band gap calculated by LDA, Eg DA \ energy band gap corrected by GWo, Ef w °; 
cohesive energy E" m , (E™) obtained with respect to nonmagnetic (magnetic) free atom energies; the C-H bond energy, E 1 q 1 1 h 
(E"c-h) obtained with respect to nonmagnetic (magnetic) free atom energies; photoelectric threshold (work function), $; 
in-plane stiffness, C and Poison ratio, v. 
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the stability of graphane structure. 

An exceptional feature of graphane is related with its 
interesting charging configuration. While graphene has a 
covalent bonding, upon the saturation of each C atom 
by a single H atom, a charge of ~ 0.1 electrons 
is transferred from H to C. At the end, the negatively 
charged bilayer of carbon atoms becomes sandwiched be- 
tween positively charged H layers. Graphane having this 
charging has the photoelectric threshold (work function) 
calculated to be $=4.97 eV, which is 0.2 eV larger than 
that of graphene. In Table |T] we presented the calcu- 
lated values related with structure, energy bands, photo- 
electric threshold etc. of graphane and graphene for the 
sake of comparison. To be complete Table [I] also includes 
calculated elastic constants such as in-plane stiffness [57] 

C [= -j ( d g ^ r ) , where Et is the total energy of the sys- 
tem under strain, e is the uniaxial strain and A is the 
area of the unit cell] and Poisson's ratio |25] v. 



IV. GRAPHANE NANORIBBONS 

The nanoribbons cut from a 2D graphane are struc- 
tures providing important features for various techno- 
logical applications. Two major families of graphane 
NRs are distinguished depending on their orientations; 
namely armchair, and zigzag graphane NRs. Apart from 
the orientation, graphane NRs are characterized by their 
widths. For armchair graphane NRs, the width, N, is de- 
fined by the number of C-C dimers in the unit cell which 
are parallel to the axis of the nanoribbon; for zigzag 
graphane NRs, N denotes the number of zigzag carbon 
chains along the nanoribbon axis. The electronic struc- 
ture changes depending on whether the dangling bonds 
of carbon atoms at the edge are free (bare) or passivated 
by single H atoms. In Fig. [2ja) and [2^b) we show the 
atomic and electronic structure of armchair and zigzag 
graphane NRs with N—ll and N=6, respectively. 

Bare armchair graphane NR is a NM indirect band 
gap semiconductor, since the electrons at the dangling 
orbitals are paired. The band gap of 2D graphane is 
reduced due to the bands of edge states of dangling bonds 
occurring in the band gap. The charge density analysis 



of these edge state bands presented in Fig. |2ja) clearly 
shows that they are localized at the edges of the NR. 
Upon hydrogen passivation of the threefold coordinated 
edge atoms, the edge state bands in the gap disappear 
and the direct band gap opens up. 

Similar to armchair graphane NRs, the dangling bonds 
of the threefold coordinated carbon atoms at both edges 
of zigzag graphane NRs give rise to edge states in the 
band gap. Owing to relatively large distance between ad- 
jacent dangling bonds and hence their relatively smaller 
coupling, they have a minute dispersion. Charge den- 
sity isosurfaces of these edge states confirm their local- 
ization at the edges. Also, owing to their weak coupling 
the spins of adjacent dangling bonds cannot be paired. 
Consequently, each carbon atom at the zigzag edge at- 
tains a magnetic moment of 1 Bohr magneton (/is). The 
ordering of these magnetic moments is however dramat- 
ically different from that of the zigzag graphene NRs. 
The significant coupling between two spin states at differ- 
ent edges gives rise to AFM ordering along the graphene 
edge. 

The physical origin of the magnetic interactions be- 
tween the edges of graphitic fragments was treated 
earlier. [31 For zigzag graphene nanoribbon, collective 
alignment of the magnetic moments along the same rib- 
bon edge through the ferromagnetic (FM) interaction, 
but AFM ordering between opposite edges is attributed 
to magnetic tail interactions. The magnetic ordering of 
a wide zigzag graphane NRs (~ 2 nm) are examined in 
a supercell comprising two unit cells, whereby antiferro- 
magnetic ordering between adjacent dangling bonds at 
the same edge is allowed. We consider four different or- 
derings, such as AFM- AFM (where, respectively, the or- 
dering of magnetic moments at the atoms located at dif- 
ferent and same edge are antiferromagnetic) , FM-AFM, 
AFM-FM and FM-FM as described in Fig. [3] Here the 
difference charge density is defined as the Ap — —pW) 
where ptw j s spin up (down) charge. Since the magnetic 
interaction between two edges vanishes, AFM-AFM and 
FM-AFM orderings have the same total energy corre- 
sponding to the ground state. These orderings are 125 
meV energetically more favorable than the AFM-FM and 
FM-FM orderings. From this analysis it is revealed that 
AFM spin alignment between adjacent atoms at the same 
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(a) 




FIG. 2: (Color online) (a) Atomic structure of bare arm- 
chair graphane NR having iV=ll. The double unit cell of 
the nanoribbon is delineated by red/dashed lines with the 
lattice constant 2a. Large/black and small/light blue balls 
indicate carbon and hydrogen atoms. Energy band structure 
corresponding to the armchair graphane NR and charge den- 
sity of selected bands are shown in the panels at the righthand 
side, (b) Atomic structure of bare zigzag graphane NR having 
iV=6 with double unit cell delineated by red/dashed lines and 
with lattice constant 2a. Energy band structure and isosur- 
face charge density of selected states corresponding to zigzag 
graphane NR are indicated. Bands shown by red/dotted lines 
are derived from edge states. Zero of energy is set to the Fermi 
level, shown by dash-dotted lines, of the nanoribbons with H- 
passivated edges. In spin polarized calculations double unit 
cell is used to allow antiferromagnetic order along the edges. 



edge is preferred by the zigzag graphane NRs. FM order- 
ing between two spins has 62.5 meV higher energy. For 
zigzag graphane NR with N=6, the magnetic coupling 
between two edges is negligible and hence either AFM- 
AFM or FM-AFM orderings have the same energy and 
they are ground states. However, when the width of the 
nanoribbon is less than 12 A the degeneracy of AFM- 
AFM and FM-AFM states is broken. The bare zigzag 



graphane NR is an indirect and antiferromagnetic semi- 
conductor with an indirect band gap relatively smaller 
than that of 2D graphane. However, it becomes NM, 
direct band gap semiconductor upon passivation of the 
dangling bonds. Also the magnetic edge states disappear 
and band gap becomes larger. 

The value of the band gap of graphane NRs can differ 
from that of parent 2D structure due to a combined ef- 
fects of quantum confinement, edge state bands as well as 
folding in the direction perpendicular to the NR axis. We 
note that the band gap of graphane NR depends on its 
width given by N. However, the quantum confinement ef- 
fect is not as emphasized as that found in graphene NR. 
In Fig. [3] H-passivated armchair and zigzag graphane 
NRs display a quantum confinement effect, namely the 
band gap reduction with increasing width or N. This 
behavior is fitted to an expression as: 

E g a P (N) = 3.42 eV + aexp(-A^) (1) 

Here a and j3 are fitting parameters. The values of a and 
(3 are found to be 1.18 eV and 0.19 for zigzag graphane 
NRs (2.15 eV and 0.14 for armchair graphane NRs). The 
band gaps of both types of graphane NRs go to 3.42 eV 
as N — > oo. Considering the GWo corrected value of the 
band gap, 5.97 eV, the above fitting can be expressed 
as Eg ap (N) = 5.97 + acxp(— N/3) assuming that the 
same scissor operation is applicable for graphane NRs. 
Apparently, the band gap values of graphane NRs cal- 
culated using DFT arc underestimated and hence their 
values are expected to occur ~ 2.5 eV larger than pre- 
sented in Fig. |4j Unfortunately, GWo self-energy cor- 
rections cannot be performed for graphane NRs due to 
large computational time. We note that the variation 
of band gap is relevant for N < 22 for zigzag graphane 
NRs and N < 30 for armchair graphane NRs. In Rcf.|20 
the band gaps of hydrogen saturated graphane NRs were 
calculated for 6 < N < 15 using Generalized Gradient 
Approximation (GGA). In the present study band gaps 
are calculated for 10 < N < 40 using LDA. The hydro- 
genated graphane NRs included in both studies, namely 
those 10 < N < 15, are in good agreement. Largest 
band gap difference of 0.1 eV found between two studies 
is attributed to different approximations (i.e LDA versus 
GGA) in representing exchange-correlation potential. 

V. FUNCTION ALIZATION OF GRAPHANE 
NRS BY ADATOMS 

Adsorption of adatoms is widely used and an effi- 
cient way to provide new functionalities to structures 
in nanoscale applications. [32 -34 Adsorption of C, O, Si, 
Ge, Pt, V, Fe and Ti adatoms on graphane nanoribbons 
is investigated by using (1x1x6) supercell of H-passivated 
armchair graphane NR having N=8. Adsorption geom- 
etry is determined by calculating the lowest binding en- 
ergy corresponding to the optimized structures. To this 
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FIG. 3: (Color online) Total energies of possible magnetic 
orderings at the edges of bare zigzag graphane nanoribbons. 
Calculations are performed in double unit cell delineated by 
red/dashed lines. Spin up and spin down states are shown by 
green/dark and grey/light isosurfaces of the difference charge 
density, Ap. 

end, the binding energies of adatoms at four initial sites 
described in Fig. |5] These are edge site E; hollow site 
HO, top site of the carbon TC and top site of the hydro- 
gen atom TH in the middle of NR. The stable binding 
sites are determined upon structure relaxation. In the 
same figure the positions of adatoms after relaxation are 
also shown. Binding energies (Ef,) of these atoms are 
calculated by using the expression 

Ejy = —Et [NR + A] + E T [NR] + E T [A] (2) 

in terms of the total energies of the optimized structure 
of graphane NR with adatom, Ex [NR + A] ; bare NR, 
Et[NR] and free adatom, Et[A}. Owing to the large 
graphane nanoribbon unit cell, k-point sampling is done 
only along the axis of the ribbon. All the energies are cal- 
culated in the same unit cell and obtained from the lowest 
ground-state total energies (either magnetic or nonmag- 
netic). Initially, at least 2 A distance between adatom 
and outermost graphane atom is provided to find the re- 
laxation site of adatom on the ribbon. 

An atom adsorbed on graphane NR may give rise to 
resonance states in the valence and conduction bands, 
and also localized states in the band gap. Owing to the 
periodic boundary conditions, the energy levels associ- 
ated with these adatoms may form energy bands. Since 
we consider a supercell consisting of six unit cell of H- 
passivated armchair graphane NR, the adatom-adatom 
distance is large and prevents coupling between adjacent 
adatoms. In this respect, the fiat bands associated with 
adatoms can mimic either resonance or localized states 
relevant for dilute doping. 

In Table [IT] we present nearest neighbor distances be- 
tween adatom and H (djy) and C (dc), binding energy 
(Eb), total magnetic moment of the system (fir), charge 
transfer between adatom and adsorbate (p*). Since the 
inversion symmetry of graphane can be broken by the 
adsorption of an atom, a net electric-dipole moment can 
be induced thereof. Thus electric dipole moment values 
calculated in the direction normal to graphane NR sur- 
face p, induced by adatom are also listed in Table |TT] 



FIG. 4: (Color online) Variation of the energy band gap of H- 
passivated zigzag and armchair NRs of graphane as a function 
of width given by N. The variation of the band gap with N 
is fitted to the curve given by continuous line. (See text) 
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FIG. 5: (Color online) Schematic representations of possi- 
ble positions of adatoms on a large H-passivated armchair 
graphane nanoribbon. Positions of adatoms obtained after 
the structure optimization through minimization of total en- 
ergy and forces exerting on the atoms are also shown. 



Photoelectric threshold of the adatom+graphane system 
perpendicular to the plane of NR, is calculated by 
the difference of the electrostatic potential at distances 
where the gradient of it is negligible and the Fermi level 
of the adsorbed system. It should be noted that both p 
and <& is relevant only for uniform coverage of graphane. 
For a dilute impurity, these calculated values converge to 
that of armchair graphane. When the width of armchair 
graphane NR is large enough, one expects the photoelec- 
tric threshold, and dipole moment, p, to converge to 
4.97 eV and el, respectively. |2T] 

We found that some atoms, such as C, O, V and Pt 
that are adsorbed on the H-passivated armchair graphane 
NR surface show a tendency to remove H atoms from 
graphane surface. Since the single atoms of C, O, V 
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TABLE II: Summary of the calculated quantities for adatoms adsorbed on a H-passivated armchair graphane NR. The first 
and second lines in each row associated with a given adatom adsorbed to edge site and the sites near the center of the graphane 
NRs, respectively. cLh, the adatom-nearest hydrogen distance; dc, the nearest adatom carbon distance; Eb, adatom binding 
energy; fir, magnetic moment per supercell; p*, excess charge on the adatom (where negative sign indicates excess electrons); 
<3>, photoelectric threshold (work function); p, dipole moment calculated in the direction normal to graphane NR surface; Ei 
energies of localized states associated with adatoms. Localized states are measured from the top of the valence bands in eV. 
The occupied ones are indicated by bold numerals and their spin alignments are denoted by either f or 4- Up to first seven 
adatom states of Ei are shown. 
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and Pt have rather high binding energy, it seems possible 
to create graphene domains on graphane by striping H 
atoms but covering the domain by adsorbed atoms. In- 
dependent of their initial position, C and O atoms are 
adsorbed to carbon atoms of graphane by replacing H 
atoms. A typical binding configuration of C and O is 
indicated as Ml in Fig. [5] Adsorption of C and O atoms 
occur at the middle of nanoribbon with strong binding 
energies of 5.89 and 3.99 eV, respectively. The resulting 
system is NM. Binding structures of Pt and V, which also 
remove H atom from the host C atom of graphane, are 
also shown in the same figure. 

Although H has very low binding energy (and hence 
is not shown in Table |n| , it may be important for strip- 
ing H atoms from graphane. Upon all geometry relax- 
ations from different initial conditions, single adsorbed 
H atom prefers to bind on top of H atom of graphane 
and forms H2 molecule which is weakly bound to the 
nearest C atom. The binding energy of H 2 molecule is 
calculated to be ~ 36 meV. 



VI. VACANCIES IN GRAPHANE NRS 

In our recent study |2T], we showed that a graphane, 
which has a NM ground state, can be made magnetic 
simply by removing hydrogen atoms from the uniform 
hydrogen layers adsorbed on its both sides. It was found 
that the magnetization depends on whether the defect 



region is one-sided or two-sided. It was also shown that 
in certain circumstances remarkably large magnetic mo- 
ments can be attained in a small domain on the sheet of 
graphane. Experimentally, removal of surface hydrogen 
atoms by using laser beam resonating with surface-atom 
bond[35|, ionic vapor application [3 6 and subnanometcr 
Pt clusters[37 are reported before. Since it is revealed 
that the double-sided vacancies give higher magnetic mo- 
ments we focus on these type of vacancies on the nanorib- 
bons. 

In this section we consider six different types of vacan- 
cies in a bare armchair graphane NR with iV=15. Upon 
desorption of a single hydrogen atom, as shown in Fig. [6] 
local bonding through sp 3 -hybrid orbital is retransformed 
into planar sp 2 - and perpendicular p z -orbitals. At the 
vacancy site one electron accommodated by the dangling 
p z -orbital becomes unpaired and hence contributes to the 
magnetization by 1 /j,b- We also calculated that a sin- 
gle H-vacancy located at the center of the NR has 70 
meV lower energy than that located near the edge of the 
ribbon. 

The situation is rather intriguing for the magnetiza- 
tion of large domains of hydrogen vacancies. As in the 
case of single H-vacancy, leading to transformation of sp 3 - 
bonding into planar sp 2 -bonding, triangular A, and rect- 
angular □ shaped large graphene islands can be created 
by removing H atoms hence by creating H-vacancy do- 
mains. As shown in Fig. [6] hydrogen desorbed triangle- 
shaped island consisting of n carbon atom at its each edge 
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is defined as An-type H- vacancy. Therefore, a An- type 
vacancy domain can be obtained by removing n(n+ l)/2 
H atom from top side and n(n — l)/2 H atom from 
down side of the graphane nanoribbon. This means that 
n(n + l)/2 p z -orbital electrons freed from A- type and 
n(n — l)/2 p z -orbital electrons freed from B-type sub- 
lattice C atoms. According to the Lieb's theorem [35] 
net magnetic moment of such a system can be calculated 
from the difference between the sublattice atoms. For 
triangle vacancies this difference is simply equal to n and 
thus the An-vacancy domain has net magnetic moment 
with n^B- However, the D-type domains does not yield 
large magnetic moments and have values between and 
2 \ib- It is seen that magnetic interactions in double- 
sided triangular and square H-vacancy domains resulting 
in net magnetic moments are relatively straightforward 
and are in good agreement with Lieb's theorem. 

We also consider CH- and C2H 2 -vacancies, the syn- 
thesis of them are relatively easier than creating only 
H- vacancies. In our study, these vacancies are first cre- 
ated in a NR by removing involved atoms and subse- 
quently their structures are optimized. Calculated val- 
ues of magnetic moments of the structure are presented 
in Fig. [6] After the creation of a CH-vacancy on a 
graphane nanoribbon, geometric structure rearranges it- 
self via J ahn- Teller- like distortion in the lattice. Charge 
density plot in Fig. [6] shows the unpaired up and down 
spin electron states located on the region surrounding the 
vacancy. It is seen that a bare armchair graphane NR, 
which is NM semiconductor when it is defect free, at- 
tains permanent magnetic moment in the presence of CH- 
vacancies. While CH-type vacancy makes small changes 
in the geometric structure, Stone- Wales type deforma- 
tion occurs after formation of the C2H2-vacancy. Since 
electrons are paired in C2H2-vacancy, the ground state is 
again NM and thus there is no visible difference charge 
density in Fig. [6] From the point view of Lieb's theo- 
rem, removal of a CH couple results in an absence in the 
number of p 2 -orbital electron belonging to a A-type (or 
B-type) sublattice. Upon the creation of C2H2-vacancy, 
both the A- and B-type sublattices lose one electron and 
thus the difference is zero resulting in a NM ground state. 
Our results are consistent with the recently reported val- 
ues for CH- vacancies in infinite sheets of graphane. |39| 



VII. EDGE ROUGHNESS 

Earlier, it was shown that the edge roughness of 
graphene NRs can affect their electronic [40 -42 and 
mechanical|43 properties. In particular, it was revealed 
that periodically repeating edge profiles can be treated 
within the superlattice structure, which result in confine- 
ment of spin states in zigzag NRs. [44 Band alignments of 
these superlattice structures have been studied both ex- 
perimentally and theoretically. [45| Here we investigated 
the effect of periodically repeating edge roughness of a 
bare zigzag graphane NR. While one edge of NR is kept 



single H W 


°0 o" C ^t =1 


|A2dcf«tfO #0 « °Ll f 2 


| A3 defect^" ° #r [1=3 












CO 9 o e O 9 


o to O O 


O O 


^* ^o o* 


D ft C. An ft O 


O o o o o o 


CO « 


o a a o 


o o o o o o o 


o *(oj .0 O 0-0 


o a o 
o o o 


• o o o 
o a o o 


O 

o o _o fo O O 


o o a o a 
o e o o o 


o a o 


• o o o 


o o o e o o o 




a o e 


■©7 Ihf J* 


O Q O O 


o o o o o o 


o o o 


o o o o 


O O O O O o 


co a o a o 


o 


■o o-a o 


O O O O O 


o o o to o o 


o 


« O 


o o o o o e o 


* O O O O O 


O O O O o • 


O O O O 


o o o a o o o 


|D 2 defect 


° ° u-2 


1 \ft A ^ 

|cndcf lc t|°. ,<3 n=i 


1 1 A ft 

iQHj defect^ ° a c [i=0 


o to 


o 


a n An a n a 


o o o o o o c 


O J O O 
CO C O 


o o o o o • o 


o o o o o o o 


o a • co o 


Co O O go 


O o c 


o * o o o « o 


O O O 3 • o O 

O O c • o 


o oo 


or o 


o « o o o o o 


a o o o o 9 


a o on a o 


o o • o • o • 


o o o ° • o 


CO 


O O 


o o o o o • a 


to to o c o 


oo oo a o 


o a • o a o • 


o o o o o * o 


o 


O O 


o p-a p-t p-4. 


o o o o o-o o 


o o o o o o 


o o • o • c o 


O O O O O o 



FIG. 6: (Color online) Atomic structures corresponding to 
single-H, double-sided triangular shaped A2 and A3, double- 
sided rectangular shaped D2, CH and C2H2 vacancies and re- 
sulting difference charge density, Ap = p"' — , surfaces for 
a bare armchair graphane NR. Large/black and small/light 
blue-gray balls indicate C and H atoms, respectively. Only 
a small part which includes vacancy region and its nearby 
atoms, of the armchair graphane NR with 7V = 15 is shown. 



flat, periodically repeating undulations are carved at the 
other edge. The same structure can be viewed as pe- 
riodically repeating heterostructures of wide (N=8) and 
narrow (N=6) segments of zigzag graphane NRs. Atomic 
and electronic band structure of this superlattice, which 
mimics the edge roughness are shown in Fig. [7] Here we 
note that infinite bare zigzag graphane NR with N=8 
has a band gap of 0.8 eV, while the band gap of infinite 
bare zigzag graphane NR with N=6 is relatively wider 
and 0.9 eV. As pointed out in Sec. IV, wide band gap 
of H-passivated graphane NRs is reduced to 0.8-0.9 eV 
because of edge states of unpassivated dangling bonds of 
bare graphane NRs, which appear in the band gap. These 
edge states in Fig. [7] of the superlattice occur as several 
closely lying flat bands below and above the Fermi level 
separated by a superlattice band gap of ~0.8 eV. The 
character of various bands are revealed by plotting the 
charge density isosurfaces of various states. Charge den- 
sity isosurfaces of the highest filled spin up and spin down 
edge state indicate the localization at the flat edge. As for 
the charge distribution of the lower lying valence states 
is rather uniform along the narrow region. The states of 
the flat band at the edge of the conduction band is con- 
fined in the wide segment of the superlattice since they 
cannot find a matching state in the narrow region. As 
a proof of the concept, it is shown that electronic states 
can be confined at specific regions of periodic edge rough- 
ness. Of course, confinement, superlattice band gap etc. 
depend on the structural parameters of superlattices and 
need detailed investigation. 
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FIG. 7: (Color online) Energy band diagram and band pro- 
jected charge density isosurfaces of various states for bare 
zigzag graphane NR including edge roughness. The band gap 
between edge states are shaded yellow/gray. Zero of the band 
energy is set to the Fermi level. 



VIII. DISCUSSIONS AND CONCLUSIONS 

Research on recently synthesized graphane revealed in- 
teresting electronic and magnetic properties of this two 
dimensional honeycomb structure. Graphane attains a 
sp 3 -like bonding through covering of its both sides by 
hydrogen atom. This bonding is rather different from 
the sp 2 -bonding of graphene and attributes a number of 
additional properties to graphane. For example, in con- 
trast to scmimctallic graphene with linear band crossing 
at the Fermi level, graphane is a nonmagnetic, wide band 
semiconductor . 

The armchair and zigzag graphane NRs of graphane 
display important number of properties and hence con- 



stitute basic structures to fabricate various devices. Both 
zigzag and armchair graphane NRs are wide band gap 
semiconductor when their edges are passivated with hy- 
drogen. The band gaps vary exponentially with their 
widths. For narrow graphane NRs the band gap is large 
due to quantum confinement effect, but approaches to the 
band gap of graphane as the width increases. Unpaired, 
unpassivated dangling bonds at the edges have 1 Bohr 
magneton magnetic moment and have antiferromagnetic 
coupling with adjacent dangling bonds. These graphane 
NRs can be functionalized to attain additional properties 
through H-passivation of their edges, adatom adsorption, 
vacancy creation, edge profiling and superiattice forma- 
tion. In particular, graphane NRs attain magnetic mo- 
ment through the creation of H vacancy at the surfaces 
of graphane. This property can be utilized to achieve in- 
teresting functionalities through decoration and pattern- 
ing of H-vacancies on the graphane NRs. The possibility 
of generation of large magnetic moments at small do- 
mains of H-vacancies makes graphane based structures 
promising for data storage and nanospintronic applica- 
tions. These functionalities can be further extended by 
adsorbing adatoms to carbon atoms deprived from hy- 
drogen. 

In conclusion, the present study demonstrates that 
graphane NRs can be an important basic nanomaterial 
and presents interesting properties for future technologi- 
cal applications. 
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